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Nous avons fait une analyse en QCD, a NLO l'approximation, de la combinaison des data de 
cross-sections profondes non-elastiques SLAC-BCDMS-NMC avec determination de a 3 . Nous 
montrons que la valeur obtenue pour a s depend de la procedure statistique du traitement des 
erreurs systematiques. L'ajustement des donnees avec la prise en compte complete des correla- 
tions point-a-point donne une valeur a s {Mz) = 0.1183 ± 0.0021(e:rp.) ± 0.0013(echellederen.) 
compatible avec les mesures de LEP et la moyenne mondiale. Nous avons extrait d'une facon 
independante des modeles les contributions "high twist" aux fonctions de structures Ft, et 
F2: la contribution de "twist-4" a Fl est en qualitatif accord avec les predictions du modele 
"renormalon" infrarouge; celle de " twist-6" a Fl tend faiblement vers les valeurs negatives et 
celle de "twist-6" a F2 vers les valeurs positives, quoiqu'elles soient toutes les deux compatibles 
avec zero compte tenu des erreurs. 
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We perform a NLO QCD analysis of the nonsingiet part of the combined SLAC-BCDMS- 
NMC data on inclusive deep inelastic cross section with the extraction of a„. We show 
that the value of ct s obtained in the analysis is sensitive to the statistical inference procedures 
dealing with systematic errors on the data. The fit with the complete account of point-to-point 
correlations of the data gives the value of a s (Mz) = 0.1183±0.0021(ea;p.) ± 0. 0013 (r en. scale) 
that is compatible with the LEP measurements and the world average. Model independent 
x-shape of high twist contributions to the structure functions Fl and F2 is extracted. Twist-4 
contribution to Fl is found to be in qualitative agreement with the predictions of infrared 
renormalon model. Twist-6 contribution to Fl exhibit weak trend to negative values, and 
twist-6 contribution to F2 - to positive values, although both are compatible with zero within 
errors. 

Data on deep inelastic scattering (DIS) of charged leptons off fixed targetsS'H are an unique 
source of information about nucleon structure and value of strong coupling constant a s . These 
data are obtained using high integral luminosity samples and their typical statistical errors 
are 0(1%). However the experimental uncertainties are dominated by systematic errors that 
typically are 2-3 times as statistical ones. Systematic errors are more difficult to be accounted 
since they are correlated and an estimators which involve correlated data are complicated. This 
is the reason, why in many analysis systematic errors are accounted using simplified approach 
and/or partially (or even completely) ignored. The aim of our present study is to perform QCD 
analysis of the data from Refsi'la with a particular attention to thorough account of point- 
to-point correlations due to systematic errors. A reliable estimate of a s implies the study of 
possible influence of high twist (HT) contribution to the scaling violation. In our analysis we 
perform simultaneous and model independent extraction of HT contribution to Fl and F2 and 
study their correlations with the a s value. 

To allow for extraction of Fl we analyzed the data on cross sections separated by the beam 
energies instead of merged data on F2. We imposed a cut of x > 0.3 to prevent additional 
uncertainties due to a poorly known gluon distribution. This cut leaves data that to a good 
approximation can be described by the pure nonsingiet structure functions, which essentially 
reduces the number of fitted parameters. The cut x < 0.75, omitting the region where the 
binding effects in deuterium are large was also imposed in the analysis. The total number of 
the data points (NDP) left after the cut is 1243; the Q range of the data is 1 GeV 2 < Q 2 < 
230 GeV 2 ; the total number of independent systematic errors is 47. 

The QCD input leading-twist (LT) structure functions of the proton and neutron were 



parametrized at the starting value of Qq = 9 GeV 2 as follows: f\ 

Ff' n (x,Q ) = A Pin x a ^(l-x) b ^, 

and then were evolved over the region of Q 2 occupied by the data in NLO QCD approximation 
in the MS factorization scheme with the help of the code used earliercl The final formula for 
the structure functions used in the fit, with account of the twist-4 contribution was chosen in 
an additive form 

F 2 ^ HT (x, Q) = F™'™%, Q) + Hjf( X )^, 

where F 2 ^'^'™ C (x, Q) are given by NLO QCD with the account of target-mass correction! 

The functions H 2 p j^\x) were parametrized in a model-independent way: their values at x = 
0.3, 0.4, 0, 5, 0.6, 0.7, 0.8 were fitted; between these points the functions were linearly interpolated. 
The data on differential cross sections were fitted using the formula 

d 2 a _ 4vra 2 (s-M 2 ) \( (Mxy) 2 \ HT ( m 2 \y 2 f^ HT ^ HT 



dxdy Q A 



where s is total c.m.s. energy, mi is scattered lepton mass and y is lepton scattering variable. 
To take into account the point-to-point correlations of the data points we minimized a 



functional 



x 2 = (fi/tx - yi)Eij(fj/£ K - yj] 



where K runs through the data subsets obtained by separation of all analyzed data on exper- 
iments and targets, and i through data points within these subsets. The matrix is the 
inverse of the covariance matrix Cy = bijOiOj + fifj(sf • s^), and each vector sf includes all 
independent systematic errors for the K-th. data subset. The other notations are: yi = the 
measurements; <7j = the statistical errors; fa = the theoretical model prediction depending on 
the fitted parameters. The normalization factors £ were fitted for old SLAC experiments and 
fixed at 1 for BCDMS, NMC, and SLAC-E-140. 

The results of this fit are given in Fig. 0. In view of large errors of H^(x) we imposed 
the constraint H^(x) = H£(x) in the final fit. The statistical quality of the fit is accept- 
able: x 2 /NDP= 1255/1243. The value of strong coupling constant obtained is a s (Mz) = 
0.1170 ± 0.0021 (stat+syst) that correspond to A^L = 337 ± 29(stat+syst) MeV or A^L = 
301 ± 30(stat+syst) MeV. Pure statistical error of a s (Mz) is 0.0011, which gives only small con- 
tribution to the total experimental error. The average bias of the the fitted function against data, 

calculated as /(/ — y)j 'y 'a 2 + f(s) 2 ^j is 0.07 that is within its possible statistical fluctuation. 

The correlation of a s with the HT contribution to F2 is very large (typical values of correlation 
coefficients is about -0.9). This means that separation of logarithmic and power effects in the 
analysis of scaling violation, which is based on the SLAC-BCDMS-NMC data without Q 2 cut, 
is unstable under various assumptions. In particular, the complete account of point-to-point 
correlations of the data leads to a shift of the a s value by about 3 standard deviations from the 
results obtained using a simplified statistical inference procedurei Since the HT contribution 
and the value of a s are strongly anticorrelated, the increase of a s is accompanied by a decrease 
of HTj] The total effect on the HT magnitude is about a factor of 3/4, as compared with the 
results of RefR At the same time Hi is almost uncorrelated with a s , i.e. its value is less model 



a We checked that extra polynomial-type factors do not improve the quality of the fits. 
6 This effect was also recently observed in the analysisp where a s (Mz) was fixed at 0.120. 
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Figure 1: The fitted values of H2(x) and Hl(x) (full circles - the fit with constraint H£(x) = H£(x), empty 
circles - the unconstrainted fit). The curves are predictions of IRR model for Q 2 — 2 GeV 2 ,n/ = 3 (full lines) 

and Q 2 = 9 GeV 2 ,n/ = 4 (dashed lines). 
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Figure 2: The dependence of twist-6 contributions to Fl and F2 on x. 



dependent. The predictions of infrared renormalon (IRR) model a are also given in Fig. g. The 
normalization factor A 2 was chosen in the universal form: A 2 = — ^^-[A^ e~ c where Cf = 4/3, 

(3) 

C = -5/3, /3 = 11 - 2/3n f , A R = A^ = 337 MeV, as obtained in our analysis. One can see 
that the model qualitatively describes the data on Hl(x) and there is evident discrepancy with 
the data on H2{x). 

We checked how much the analyzed data are sensitive to the twist-6 contribution to Fl 
and F 2 . For this purpose we added to F 2 HT or F^ T the terms H^ 2 (x) 1 q1 v , where functions 

Hl\(x) were the same for proton and deuterium, and parametrized similarly to Hi j2 (x). The 

fitted values of H^ 2 (x) are given in Fig. [| One can observe the trend to the negative values 



at highest x for H L i \x) and trend to the positive values for H^'(x). However the statistical 
significance of the deviation of H^\(x) from zero is not very large. In addition, the correlation 
of H^\x) with a s and H 2 (x) is extremely strong. Summarizing these observations, we can 
conclude that the observed deviation of H L i \x) and H^\x) off zero can be considered as 
qualitative only. 

Any renormalization group analysis, which is based on the finite number of perturbative 
series terms is sensitive to the choice of renormalization scale. This dependence can be used for 
the rough estimate of the higher orders terms effect that are not accounted for in the analysis; 
if one includes all terms this dependence should vanish. For estimate of renormalization scale 
uncertainty in DIS analysis this scale usually is being chosen as KrQ 2 and the value of Kr 
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Figure 3: The dependence of a a (Mz) on the renormalization scale. Open circles correspond to the fits with the 
HT contributions fixed, full circles - to the fits with the HT contributions released. 




Figure 4: The values of HT contributions for Kb. = 0.25 (open circles) and for Kr = 1 (full circles). 



is varying in the range of 0.1-4li<J We studied the dependence of a s on the value of Kr in two 
cases: for the HT contributions released in the fit and for the HT contributions fixed at the values 
obtained in the fit with Kr = 1. The results are given in Fig. ||. One can see that the dependence 
of a s on the renormalization scale is different in the two cases. The reason of this difference 
is that the H^' (x) are strongly correlated with a s ; when one changes Kr, the changes in the 
QCD evolution kernel can be absorbed into the additional power-like contribution. This effect, 
illustrated in Fig. |4], can be considered as indirect indication that the fitted values of H^ix) can 
include not only genuine power corrections, but also effectively account for higher orders QCD 
terms. Any way the the value of a s is quite stable against renormalization scale variation, if HT 
are released. At the same time the a s errors are significantly larger for this case, which is also 
a consequence of the large correlation of a s with HT, i.e. one can say that the uncertainty in 
the ot s value due to the renormalization scale choice is partially hidden in the total experimental 
error. The spread of the a s value in the fits with Kr =0.25-4 and HT released is 0.0026. With 
the account of this spread we obtain a s (Mz) = 0.1183 ± 0.0021 (exp.)±0.0013(ren. scale). 
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